Weak disorder: anomalous transport and diffusion are normal yet again 
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We have carried out a detailed study of the motion of particles driven by a constant external 
force over a landscape consisting of a periodic potential corrugated by a small amount of spatial 
disorder. We observe anomalous behavior in the form of subdiffusion and superdiffusion and even 
subtransport over very long time scales. Recent studies of transport over slightly random landscapes 
have focused only on parameters leading to normal behavior, and while enhanced diffusion has been 
identified when the external force approaches the critical value associated with the transition from 
locked to running solutions, the regime of anomalous behavior had not been recognized. We provide 
a qualitative explanation for the origin of these anomalies. 
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Solid state surfaces frequently present periodic poten- 
tials marred by some disorder. Herein we show that an 
overdamped particle moving over such a potential in one 
dimension (ID) may exhibit anomalous behavior in the 
form of superdiffusion, subdiffusion, and even subtrans- 
port. Although we cannot prove that these are steady 
state regimes, our numerical simulation data show them 
to be present over time spans of several orders of magni- 
tude. 

That diffusion of particles over both periodic and ran- 
dom surfaces lead to some forms of anomalous behav- 
ior is of course well known and continues to attract a 
great deal of attention both theoretically and experimen- 
tally [lrlTfjj. In periodic potentials with low friction, ex- 
tremely long (in time) dispersionless transport regimes 
can be observed when forces exceed a critical force @. 
Moreover, in these same systems, in both overdamped 
and underdamped regimes, the diffusion coefficient ver- 
sus the applied force presents a pronounced peak around 
the critical force that allows the coexistence of locked 
and running states [l|, 0-0] . The The enhancement is 
quantitatively larger than the free particle diffusion coef- 
ficient. This behavior has been observed experimentally 
when tracking the motion of colloidal spheres through a 
periodic potential created with optical vortex traps Q . 

The enhancement of the diffusion coefficient is even 
more pronounced when disorder is also present @. This 
phenomenon has been tested by numerical simulations on 
a surface in which a small amount of spatial disorder in 
the form of a random potential is added to the periodic 
potential [Toj . Dramatic diffusive enhancement occurs 
even for very small amounts of disorder, e.g., when the 
amplitude of the random contribution of the potential is 
as small as ~ 5% of that of the periodic contribution. 

Although dramatic, diffusive enhancement turns out to 



be only a limited aspect of the story because it is not the 
only manifestation of disorder. Here we present a range 
of additional anomalous transport and diffusion phenom- 
ena arising from weak disorder that have not been pre- 
viously noted. Our model and the behaviors it exhibits 
are inspired by [t| [l(| . 

We consider the overdamped motion of identical non- 
interacting Brownian particles moving in a ID poten- 
tial landscape U (x) following the Langevin equation 
7x(i) = —U'{x) + F + £(t). Here x is the position of 
the particle, t denotes time, 7 is the dissipation param- 
eter, F is the applied force, and £(t) is Gaussian ther- 
mal noise at temperature T . The correlation function 
of the noise obeys the fluctuation-dissipation relation 
(£(*)£(£')) = 2~/k B T5(t - t'). The potential U(x) con- 
sists of a periodic part, V p (x) = Vo cos (2irx/\ p ), and a 
Gaussian spatially random contribution V r {x) with cor- 
relation function 

T/2 / o 2 2\ 

g r {x) ee (V r (x)V r (0)) = -^cxp ^ —J . (1) 

We need to choose parameters to capture the relative con- 
tribution of each, as well as their relative amplitudes and 
length scales. The relative contributions are determined 
by the parameter a in the combination 

U(x) = (l-a)V p (x)+aV r {x), <jG[0,1]. (2) 

Furthermore, we have chosen the potential correlations 
g r (x) and g p {x) = (V 2 /2) cos (2ttx/X p ) to be equal at x = 
0, ffp(0) = ffr(0). This ensures that the total potential 
amplitude is of order Vo independently of a. Also, with 
the particular choice l r — X p the two correlation functions 
are identical up to second order in a Taylor expansion. 

The equation of motion can be rescaled into dimension- 
less form in terms of the spatial variable z = 2ttx/X p and 
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the temporal variable r = [(2ir) 2 Vo/jXp]t. This yields 

z = (1 - <r)/ p (z) + (a/X)f r (z/X) +T + n(r), (3) 

where f p and / r are the dimensionless forces arising from 
V p and V r , respectively, and r\ is the dimensionless noise. 
The dimensionless parameters are 
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Throughout this work we set T — 0.01 and X p = 2tt. 
Variations in T do not lead to any additional phe- 
nomenology. The specific choice of X p is only impor- 
tant for passage between dimensionless and dimensioned 
units. An important note is that in these variables a de- 
crease of A even for fixed a leads to an increase in the 
relative contribution of the random force. 
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FIG. 1. (Color online) Upper panel: Barriers h of one realiza- 
tion of the total potential in a spatial domain of ~ 32 periods 
with A = 0.5 and T — 1. Middle and lower panels: Numer- 
ical data for 50 different disordered potentials and N = 100 
particles in each potential. Left panels: v{r). Right panels: 
D(t). Middle panels: A = 1. Lower panels: A = 2n (regime 
of Ref. [U}). a = 0.05, r = 10 4 . 

We have carried out numerical simulations of Eq. ^ 
over a large number (100) of particle trajectories and 
a large number of realizations of the random potential 
contribution (typically 50-100) in order to calculate the 
velocity and diffusion coefficient using the prescriptions 
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The brackets (• ■ • ) denote averages over many trajecto- 
ries and potentials. The numerical procedures are en- 
tirely standard. The usual assumption is that stationary 
values are reached at long times, an assumption that we 
show here to be at the very least questionable. Repre- 
sentative outcomes of this procedure are shown in the 



middle and lower panels of Fig. [T] These outcomes qual- 
itatively capture an unanticipated range of behaviors in 
the right middle panel. 

The velocity reaches a well-defined stationary value for 
any value of the force in the figure (left panels), even in 
the regime of sharp increase around J- = 1. Also, as 
noted earlier, for the parameters used in Ref. [10|, we see 
the enhancement of the (well defined) diffusion coefficient 
around the critical deterministic force J- c , as discussed in 
that work (lower right panel). However, entirely differ- 
ent behavior is now seen in the middle right panel of the 
figure, which shows huge variations in the diffusion coef- 
ficient. Note that these variations extend over orders of 
magnitude. The question then is - what leads to these 
variations and why were they not identified in Ref. [Toj? 

The answer lies in the choice of the correlation length 
of the random potential. In fiol ] the correlation length is 
much greater than a single period of the periodic poten- 
tial, that is, the randomness is very smooth. However, 
entirely different behavior is seen when the random po- 
tential is more corrugated, as it is in our case. We see 
that the traditional diffusion coefficient is no longer well- 
defined in this regime but instead presents very strong 
fluctuations around its peak value, making a precise esti- 
mation of D questionable. Note also the very large scale 
differences of the middle and lower right panels of the fig- 
ure. These are signatures of anomalous behavior. While 
it is not clear whether a well-defined value of D would be 
obtained at much longer times beyond our computational 
reach, we have repeated these calculations for many more 
particles and potential realizations and continue to find 
this variability. 

To understand the origin of these anomalies, in the 
upper panel of Fig. [1] we plot the barriers h of the total 
potential U (z)—Tz for a particular realization of the ran- 
dom potential. The heights and locations of the barriers 
are random, and most of them exceed the thermal energy 
(dashed line). Smaller values of A lead to a greater num- 
ber and height of the barriers. As a particle moves along 
such a landscape, it must overcome these barriers, some 
of which are extremely high. Indeed, V r {z) allows for 
barriers of any height, limited only by finite system sizes 
and simulation times. Thus, even with a small amount 
of disorder the particle motion is dominated by random 
waiting times due to the dispersion of the barrier heights, 
and a few long waiting times that greatly influence the 
outcome. We go on to show that the diffusion anomaly in 
Fig. [T]is a consequence of such landscapes and that it is 
qualitatively different from a simple large enhancement 
of the diffusion coefficient. In fact, strong fluctuations in 
the usual ensemble calculation of the diffusion coefficient 
point to the fact that our system is exhibiting behavior 
reminiscent of aging or of weak ergodicity breaking [Tll - 
[l3[ over the time scales of our simulations. We approach 
the problem with this observation in mind. 

Our numerical results are collected as follows. Parti- 
cles are initially located at random positions uniformly 
extended over a region of about 1000 sites, and for each 
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we observe the times t p that it takes a single particle to 
cover the underlying spatial period X p over the course of 
its trajectory over a long time. We collect these statistics 
for many particles and many realizations of the random 
potential. If the motion of the particles is "normal" then 
following the reasoning in [Hill and also used in [13] 
leads to the average velocity and diffusion coefficient ex- 
pressions 
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Here (t p ) is the average of the t p over all trajectories, 
particles, and potential realizations and (A 2 i p ) is their 
dispersion about the average. It is of course evident that 
these expressions are only well defined if the first and 
second moments of the distribution of the t p are finite. 
We present numerical evidence that indicates that in the 
presence of a small amount of disorder in the potential 
the distribution P{t p ) of the time t p to cover a single 
period can have a power-law tail for long stretches of 
time, 



P(t p ) ~ t~ 



(7) 



This behavior leads to the observed anomalies if (3 is suf- 
ficiently small. For j3 > 3, the first and second moments 
are finite, so transport and diffusion are normal. In the 
range 2 < (3 < 3, the first moment is finite but the second 
moment diverges. This leads to a finite average velocity 
v, but the diffusion coefficient as defined in Eq. ^ di- 
verges (superdiffusion) . In the interval 3/2 < (3 < 2 
the first moment diverges and the average velocity thus 
vanishes (subtransport). The second moment again di- 
verges (superdiffusion). For 1 < (3 < 3/2 we again 
have substransport, and the diffusion coefficient decays 
to zero (subdiffusion) . In this case the particle remains 
extremely localized around its point of origin. Interest- 
ingly, we observe all of these behaviors. 

To obtain values of (3, we focus on the tails of the dis- 
tribution. If they follow a power law, exponents (3 are 
estimated, as depicted in Fig. [3J In the upper panel of 
Table U we present a set of typical results of simulations 
leading to different values of the exponent /?, obtained 
for two levels of disorder and for different forces. We ob- 
serve that, in many cases, a value (3 < 3 associated with 
anomalous behavior is found. Table Q] also shows that 
a lower level of disorder (~ 2%) leads to higher expo- 
nents. These estimated values are of course informative 
but should not be interpreted too narrowly because they 
are based on sparse statistics. For the sake of compari- 
son with Ref. |10[, we also show the case A = 2n for a 
5% level of disorder. We observe what we expected from 
Eq. §3§ when A increases, that is, the effect of the ran- 
dom part of the potential is greatly muted. This implies 
that most of the exponents are in the range f3 > 3 (nor- 
mal behavior) for the range of forces selected, which ex- 
plains why the anomalies explored herein were not found 
in Ref. Much smaller forces would need to be ex- 

plored to find robust anomalous behavior there. For our 




FIG. 2. (Color online) Log-log plot of the time distribution 
for different values of A. (a) A = 0.1 (/3 = 1.4), (b) A = 0.4 
(/3 = 1.8), (c) A = 0.8 [p = 2.6), and (d) A = 2tt (/3 = 
6.9). The parameter values are T = 1 and a = 0.05. Inset: 
Log plot of P{t p ) for cases (c) and (d). Also shown is the 
exponential distribution for a purely periodic potential with 
the same force (green online). Time r > 10 4 as needed to 
obtain reliable histograms. 



parameter choices, using A as a control parameter (keep- 
ing T fixed) also leads to f3 values covering all the possible 
regimes (see bottom panel of Table [TJ> . 
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TABLE I. Upper panel: Exponents f3 for different values of 
the force and two values of a and A. Bottom panel: Exponents 
/3 for different values of A, with T — 1 and a — 0.05 



As the value of /? increases with increasing A, the shape 
of the distribution P(t p ), particularly its tail, changes 
from a power law form to the more typical exponential 
associated with normal behavior. This occurs because 
changing A induces a change of the effective contribution 
of the random part of the potential, as already noted 
earlier. We pointed out that increasing A diminishes the 
effects of the disorder. In fact, case (d) of Fig. [2] is no 
longer a power law, as shown in the inset of the figure, 
but is instead better described as an exponential, char- 
acteristic of normal diffusive behavior. 
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FIG. 3. (Color online) Log-log plot of the time evolution of 
the velocity (left) and of the diffusion coefficient (right) for 
different values of A: A = 0.1 (solid line), A = 0.4 (dashed 
line), A = 0.8 (dotted line) and \ — 2n (dotted-dashed line). 
Parameter values: a — 0.05, T = 1. Number of potential 
realizations: 100, except for A = 27r where 20 realizations are 
sufficient. 

Finally, temporal evolutions of the velocity and the dif- 
fusion coefficient for the four different regimes of behavior 
are shown in Fig. [3] They are in qualitative agreement 
with this analysis, at least during a transient of several 
decades. The subtransport regimes of the two cases with 
the lowest values of A are clearly seen, as are the subdif- 
fusive and superdiffusive behaviors of the diffusion coef- 
ficient. 

In conclusion, we have carried out a detailed numeri- 
cal study of the motion of particles driven by a constant 
external force over a one-dimensional landscape consist- 
ing of a periodic potential modified by a small amount 



of spatial disorder. We have identified a set of dramatic 
anomalous behaviors as diverse as subtransport, subd- 
iffusion, and superdiffusion on the same surfaces as the 
driving force or the random corrugation length is varied. 
These behaviors are observed over very long time scales. 
Their asymptotic persistence behavior is not known. Ear- 
lier studies have focused only on parameters of normal 
behavior [ij [l(|, and while they have identified the oc- 
currence of enhanced diffusion when the external force 
approaches the critical value associated with the tran- 
sition from locked to running solutions, they have not 
recognized the regime of anomalous behavior. 

The regimes that exhibit these anomalous behaviors 
are identified by the correlation length l r of the random 
portion of the potential. We find anomalous behavior 
when l r < A p , the period of the periodic portion, that 
is, when the periodic potential is slightly corrugated over 
short distances. Earlier studies had focused on the regime 

= 27rA p , that is, on regimes of very long smooth vari- 
ation of the random contribution to the potential. Ex- 
periments with more corrugated surfaces than have been 
used so far [9j are clearly desirable to see the effects that 
we have identified. 

It would of course also be desirable to extend these 
studies to higher dimensions and to find an analytic char- 
acterization of these systems, one which would allow in- 
sight into the asymptotic behavior. Such an analysis has 
recently been presented for the case of a piecewise linear 
random potential [lij], but seems not yet to be available 
for the more realistic potentials considered here. 
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